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Abstract. 

We present a theoretical study of the Josephson couphng of two superconductors 
which are connected through a diffusive contact consisting of noncolhnear 
ferromagnetic domains. The leads are conventional s-wave superconductors with a 
phase difference of ip. First, we consider a contact with two domains with magnetization 
vectors misoriented by an angle 9. Using the quantum circuit theory, we find that in 
addition to the charge supercurrent, which shows a — tt transition relative to the angle 
^, a spin supercurrent with a spin polarization normal to the magnetization vectors 
flows between the domains. While the charge supercurrent, is odd in if and even in 
0, the spin supercurrent is even in ip and odd in 6. Furthermore, with asymmetric 
insulating barriers at the interfaces of the junction, the system may experience an 
antiferromagnetic-ferromagnetic phase transition for (p = n. Secondly, we discuss 
the spin supercurrent in an extended magnetic texture with multiple domainwalls. 
We find the position-dependent spin supercurrent. While the direction of the spin 
supercurrent is always perpendicular to the plane of the magnetization vectors, the 
magnitude of the spin supercurrent strongly depends on the phase difference between 
the superconductors and the number of domain walls. In particular, our results reveal 
a high sensitivity of the spin- and charge-transport in the junction to the number 
of domain walls in the ferromagnet. We show that superconductivity in coexistence 
with non-coUinear magnetism, can be used in a Josephson nanodevice to create a 
controllable spin supercurrent acting as a spin transfer torque on a system. Our results 
demonstrate the possibility to couple the superconducting phase to the magnetization 
dynamics and, hence, constitutes a quantum interface e.g. between the magnetization 
and a superconducting qubit. 
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Figure 1. (a) (Color online) S1F1F2S2 junction, where, each ferromagnet is presented 
by a node. 9 is the angle between the magnetization direction of the Fi and F2. (b) 
The S1F1DWF2S2 junction with Neel domain wall (DW), a is the angle which the 
local exchange field in the domain wall makes with the y-direction. 



1. Introduction 



The Josephson effect refers to the coherent transfer of Cooper pairs between two weakly 
coupled superconductors |1]. In a Josephson contact with a normal metal between 
the superconductors, the underlying microscopic mechanism is Andreev scattering |2] 
at the two normal-metal-superconductor interfaces which converts electron and hole 
excitations of opposite spin directions into each other by creating a Cooper pair. 
The resulting dissipationless electrical current is driven by the difference between the 
phases of the superconducting order parameters across the contact. From the fact 
that superconductivity is a coherent state of spontaneously broken U(l)-symmetry, it 
follows that the Josephson effect is a response of this coherent state to an inhomogeneity 
over the junction which is produced by the variation of the phase. An analogous non- 
superconducting effect is predicted to exist in the magnetic tunnel barrier between two 
ferromagnets with the S0(3) symmetry breaking coherent states |3l H]. In this case the 
misorientation angle 6 of the two magnetization vectors is the driving potential for a 
dissipationless spin supercurrent, similar to the exchange interaction. 

In this Article we develop the circuit theory of the superconducting spin Josephson 
effect in an inhomogeneous ferromagnetic (F) contact between two conventional 
superconductors. We show that when the F contact consists of two domains whose 
magnetization vectors enclose an angle 6, in addition to the charge supercurrent, 
a spin supercurrent will also appear. This spin supercurrent is created by the 
simultaneous existence of the superconducting states and a noncollinear orientation 
6 of the magnetization vectors. Interestingly, we find that the spin supercurrent and 
the corresponding spin transfer torque is directed perpendicular to the plane of the two 
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magnetization vectors, such that it would lead to a precession of the magnetizations 
around each other. In addition to extensive theoretical O [6l [71 E] and experimental 
P [ini EH [El Us] studies of the spin-transfer torque in F spin-valve and domain 
structures, there have been studies devoted to the spin-transfer torque in structures 
with superconducting parts [Tl[T5l[l6l[2a[I71[l8l[T9l[2nilSl[22l[23l[23 The 
essential effect is the production of long-ranged spin-triplet superconducting correlations 
by the interplay between the induced spin-singlet correlations and the noncoUinearity 
of the magnetization profile in F contact |271 [281 EHl [3Ql [SH [321 ES]. Compared to 
previous studies, we present a quantum circuit theory calculation which takes the spatial 
variation of Green's functions as well as the nonlinearity of the proximity effect fully 
into account. By this method we specifically are able to obtain the inhomogeneity of 
the spin supercurrent and to define a spin-transfer torque in non-collinear ferromagnetic 
Josephson contacts. 

Several experimental works on this triplet proximity effect have been done [311 [25] . 
Recently, Khaire et al. [36] reported the observation of the long-range supercurrent 
in Josephson junctions which is controllable by varying the thickness of one of the 
ferromagnetic domains. Also, Robinson et al. [371 EH] detected the flow of a long-range 
supercurrent in the ferromagnetic Josephson junction with a magnetic Ho-Co-Ho trilayer 
and found an enhancement of the critical currents in the antiparallel configuration of 
the Junctions with a trilayer Fe/Cr/Fe barrier. 

In analogy to the conventional charge Josephson effect, the spin Josephson effect 
has the tendency to remove the inhomogeneity of the order parameter vector of 
the spin-triplet superconducting state in the F-contact. This is analogous to the 
spin Josephson effect in contacts between two unconventional triplet superconductors, 
where the Cooper pair spin current appears in conjunction with the usual charge 
supercurrent [an SQl SB [12] . 

The spin-dependent circuit theory has already been used in Ref. [32] to study the 
density of states, and the Josephson supercurrent in S/F/S heterostructures, which are 
shown to be dependent on the configuration of the magnetization in F. Here, we further 
study the spin supercurrent and spin transfer torque in such Josephson junctions. We 
demonstrate the dependence of the charge and spin supercurrent on the phase difference 
(/? and the angle between the magnetizations 9: the spin supercurrent is an even function 
of ip and an odd function of 9, the charge supercurrent satisfies the inverse relations 
relative to the (p and 9. Further, we study the equilibrium configuration of the exchange 
field vectors as a function of the phase difference and the temperature. We obtain 
phase diagrams which show the antiferromagnetic-ferromagnetic phase transitions in 
the system which hasn't been announced in any other literatures. 

We also discuss the generality of this effect for other ferromagnetic contacts with a 
more complex inhomogeneity of the direction of the magnetization vector. In particular, 
when the ferromagnetic contact consists of an in-plane rotating magnetization vector 
between two homogeneous domains with antiparallel magnetization, we find that the 
spin supercurrent is highly sensitive to the value of the wave vector. We show that one 
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can tune the spin supercurrent acting as a spin transfer torque, by changing the phase 
difference between the superconductors or with variations of the wave vector. Also, we 
investigate the position dependence of the spin current in S1F1DWF2S2 which, to the 
best of our knowledge, hasn't been studied in any other texts. We find that the behavior 
of spin supercurrent relative to the position strongly depends on the phase difference 
between the superconductors. We extend the quantum circuit theory by the description, 
how spin-transfer torques can be calculated within the method. 

2. Model and basic equation 

We describe the basic theory and the model first for a two-domain ferromagnetic contact 
between two conventional superconductors, as is shown schematically in Fig. [1^. The 
exchange field of one domain Fi makes an angle 6 with that of the other domain F2. 
We restrict our study to the time-independent case and do not consider changes in the 
magnetic structure in this article. The generalization to the structure with a continuous 
magnetization texture in Fig. [Tb is straightforward and described in the end of this 
section. To proceed our work, we make use of the quantum circuit theory which is a 
finite-element technique for calculating the quasiclassical Green's functions in diffusive 
nanostructures [l3l HU HSl HHl HTJ HH]. In this technique, we represent each F domain 
and S reservoirs by a single node, which is characterized by an energy-dependent 4x4- 
matrix Green's function Gi, in Nambu and spin (='|", J,) spaces |17l HE]. Furthermore, 
the two nodes in the F domains are assumed to be weakly coupled to each other by 
means of a tunneling contact. In terms of its spin-space matrix components g and /, 
the matrix Green's function is written as 



We consider the equilibrium condition where a misorientation angle 6 and phase 
difference cp may drive Josephson spin and charge currents between two adjacent nodes. 
These equilibrium Josephson currents can be extracted from the matrix current defined 

as 



where gij is the tunneling conductance of the contact between two nodes and i and j 
denote the connected nodes. This approach works also, if we divide the ferromagnetic 
region into n nodes, as is necessary in the case of Fig. [H Then we have to take 
{n — l)/gij = {1/ gFiF2) ~ iXl QsiFi + ^l gs2F2)- Where, the conducting part of F domains 
is discretized into n nodes. For a S1F1F2S2 structure, n is two. The conductance of 
the tunnel barrier between Sl(2) and Fl(2) is denoted by gsi{2)Fi{2) and, gFiF2 is the 
conductance of the whole F1F2 contact. 

The matrix current obeys the following law of current conservation in matrix form 




(1) 



(2) 



(3) 



j=i±i 
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Here I^i = -Gg^uj/Si) 



f^jGi , with 6i being the electronic level spacing of the 
node, is the matrix of the leakage current which takes into account the dephasing 
of electrons and holes due to their finite dwell time in the node i, and Igi = 
i{GQ/6i) {h^CFi + /if (3-2 + hlf^a^),Gi is the corresponding matrix current representing 

the leakage caused by the spin-splitting due to an exchange field hi {Gq = e^/{2Tih) 
is the quantum of the conductance). The third term represents the matrix currents 
from the neighboring nodes i-l,i+l. From Eq. ([3]) we can find the spin-torque, e.g. in 
z-direction as r^i = Izi,i+i + Izi,i-i- 

Equation ([3]) is given for all nodes and is supplemented by boundary conditions, 
which are the values of G in the S reservoirs. We neglect the inverse proximity effect in 
the reservoirs and set the matrix Green's function in SI and S2 to the bulk values: 

6m - ^fi^ (4) 



^2 + |A'' 



where 





Ai,2= xt . '^3= , , . (5) 



Here, Ai_2 = | A| exp (±i(/?/2)(Ti are, respectively, the superconducting order parameter 
matrix in SI and S2 ((Xj denote the Pauli matrices in spin space) and u = 
7rT(2m + 1), with m being an integer, is the Matsubara frequency. The temperature 
dependence of the amplitude of the order parameter is well approximated by |A| = 



1.76Tctanh(1.74-y/ Tc/T — 1). We note that the matrix Green's function satisfies the 
normalization condition G^ = 1. 

We have solved these equations numerically by an iteration method. In our 
calculation we start by choosing a trial form of the matrix Green's functions of the 
nodes, for a given 0, T, and the Matsubara frequency m=0. Then, using Eq. ([3]) 
and the boundary conditions iteratively, we refine the initial values until the Green's 
functions are calculated in each of two nodes with the desired accuracy. Note that in 
general for any phase difference 0, the resulting Green's functions vary from one node to 
another, simulating the spatial variation along the F contact (see Ref. |19]). From the 
resulting Green's functions and Eq. ([2]) we find the matrix currents. Then, we calculate 
the charge supercurrent I and the components of the spin supercurrent vector I from 
the relations 

lij = tra^Iij , Izij = trf^a^Iij , Ix(y)i,j = ^^(^i{2)hj ? (6) 

in which tr . . . = (i7rT/2e) J2ui Tr . . . with Tr denoting the trace in Nambu-spin spaces. In 
the next steps we change to the next Matsubara frequency and use the results from the 
previous one cj^-i as initial guess. We find the respective contribution to the spectral 
currents and continue to higher frequencies until the required precision of the summation 
over m is achieved. 

In the following we scale the length of the system, L, in units of the diffusive 
superconducting coherence length ^, and use the dimensionless parameters of h/Tc and 
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t = T/Tc, as measures of the amplitude of the exchange field h and the temperature T.To 
describe a continuous domain wall, we use as parameter the wave vector Q associated 
with one full winding of the magnetization by 27r. Hence the total number of windings 
of the magnetization is given by QLp/2'K, where Lp is the length of the inhomogeneous 
region (see Fig. [1]). The currents are expressed in units of /q, the amplitude of the 
critical Josephson charge current at /i = and t = 0. 

3. Results and discussions 

From the numerical calculations we obtain the spin supercurrent with a polarization 
directed normal to the plane of the magnetization vectors. Further, we find a transition 
of the favorable configuration of the domain, from antiparallel to the parallel as the 
exchange field of the asymmetric domains increases. Also, we show that in a system 
with more complex configuration of the direction of the magnetization, the profile and 
penetration depth of the spin supercurrent are highly dependent on the number of the 
rotations that the magnetization vector has undergone across the domain wall. 

3.1. S1F1F2S2 junction 

Using the method described in Sec. 2 we have calculated the spin and charge 
supercurrents for the two domains F contact of Fig. [TK, when the exchange field 
vectors are taken to be in the x — y plane. We found that the spin supercurrent has a 
polarization which is aligned along the z axis, namely perpendicular to the plane of the 
exchange fields of Fi and F2. Our results for the dependence of the spin 1^ and charge / 
supercurrents on the misorientation angle 6, the phase difference ip and the temperature 
t are shown in Fig. [2l We found that, in general, the spin supercurrent obeys the 
symmetry relations Iz{f) = Iz{—v) and IziO) = which are the analogs of the 

relations I{ip) = —I{—ip) and I{9) = I{—9) for the charge supercurrent (see Fig. [2]a,b). 
These behaviors suggest that one can change the direction of the spin supercurrent, 
which is proportional to the induced spin transfer torque [7j, by changing the phase 
difference between two superconductors when 6 is fixed. We note that a nonzero spin 
supercurrent is provided by a noncollinear orientation of the exchange field vectors and 
the existence of the superconductivity (|A| 7^ 0), even for (p = 0. 

We define the critical spin supercurrent, I^cri'^), as the maximum of the absolute 
value of the spin supercurrent as a function of 6 for a given ip, in similarity to the 
definition of the charge critical supercurrent. We may also use a distinct definition, 
which we denote by lzmax{0), as the absolute value of the spin supercurrent for a value of 
(p which maximizes the charge supercurrent as a function of p, for a given 6. Figure. [Jt 
shows the behavior of Izcr as a function of p> for different temperatures. At a given 
temperature t, Izcr{p) shows a change of sign at a phase difference which depends on t. 
This change of sign may be recognized as the signature of a transition between and 
TT spin Josephson couplings, in analogy to the charge — vr transition in F Josephson 




Figure 2. (a) (Color online) Plot of spin supercurrent versus 9 for different values of 
tf when L/^ = 1.0, h/T^ = 5.0, t = 0.5. (inset) 9 dependence of charge supercurrent 
for the same system, (b) Plot of / versus for different values of 9 for the previous 
system, (inset) Lp dependence of spin supercurrent for the same system, (c) Critical 
spin supercurrent versus for different temperatures, which shows the appearance 
of — TT transition relative to the 0. (inset) Critical charge supercurrent versus 9. 
(d) Light dashed, dotted, dashed-dotted lines are, respectively, the maximum spin 
supercurrent when ip changes between and tt versus 0, when t = 0.1,0.5, and 0.9. 
Dark dashed, dotted, dashed-dotted lines are, respectively, spin supercurrent for the 
Lp which maximize the charge supercurrent for different temperatures for the same 
situation. 

junctions [501 ED ESj [531 EH ESI EEl EZl EH EH] . The corresponding critical charge current 
shows a — vr transition with varying 9. Note that both Izcriv) and Icr{d) have a nonzero 
value at the transition point at low temperatures, as the signature of nonzero second 
harmonic in the current-phase and the current-angle relations [60 1 [6T 1 162] . Fig. |2]i we 
have also plotted Izmax{9), which shows a change of sign at ^ = vr for all temperatures 

m- 

We have also studied the dependence of the spin supercurrent on the absolute 
value of the exchange field. The results are shown in Fig. |3l in which I^cr is plotted 
as a function of h/Tc for different (p and 6. These results show that the sign and the 
amplitude of the spin supercurrent can be also modulated by varying h/Tc, which can 
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Figure 3. (Color onliiie)(a) Iz versus h/Tc for different values of 9 when (p — tt/2, 
L/^ = 1.0, and, i = 0.1. (b) 1^ versus h/Tc for the same system but for different 
if when 9 = Tr/2. (a) and (b) are logarithmic plots which show the — tt transition 
relative to the exchange field. 



be used for further tuning of the corresponding spin transfer torque. We note that for 
strong ferromagnets with h ^ A, the spin supercurrent vanishes. This is due to the 
suppression of the amphtude of the Andreev reflection at SiFi and S2F2 interfaces in 
this hmit, which suppresses the proximity effect. 

It is also interesting to study the equilibrium configuration of the exchange field 
vectors as a function of the phase difference and the temperature. The equilibrium angle 
can be obtained by minimizing the free energy, F, of the contact as a function of 6. We 
have calculated the ^-dependence of F by integrating the spin supercurrent over 6 and 
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Figure 4. (Color online) Phase diagram of the transition of the minimum of free 
energy from 9 = tt to 9 = for asymmetric systems with gsiFi = 0-^gs2F2 and 
different lengths. The inset shows the F — 9 plane of the 3-dimensional energy plot 
when L/^ — 1 and t — 0.4. Here, _fo is the free energy of the system when the 
magnetization is coUinear. 

charge supercurrent over 

F{y,,6) = r i{ip' ,Q)dip' + f h{ip,e')de'. (7) 

JO Jo 

Our calculation shows that the exchange field vectors favor either parallel [9 = 0) 
or antiparallel {6 = it) configurations, depending on (p, t and h/Tc. The behavior of 
this superconductivity-induced exchange coupling differs for the two cases of a contact 
with symmetric barriers with Qsifi = gs2F2 and an asymmetric contact with a very 
different gsiFi and gs2F2- For a symmetric system, we have found that the coupling is 
antiferromagnetic {9 = tt) for ip = 0, but becomes ferromagnetic {6 = 0) for = vr. 
This behavior is found to hold irrespective of the values of L/^, h/Tc. 

However, for an asymmetric system it is possible to change the coupling from 
ferromagnetic to antiferromagnetic and vice versa by varying L/^, h/Tc, or t, for the 
phase difference cp = ir. In Fig. HI we have shown the ferromagnetic-antiferromagnetic 
coupling phase diagram of the system in the plane of h/T^ and t, when p = tt and for 
some different values of L/S,- This phase diagram is similar to the — vr Josephson 
couplings phase diagram of a homogenous F contact between two superconductors, see 
Ref. jH]. As we show in the inset of Fig. HI the minimum of F as a function of 9 
shifts from = tt for low values of /i/Tc to = at higher h/T^s, when L/C, = 1 
and t = 0.4. The temperature induced transition between the ferromagnetic and 
the antiferromagnetic phases is also possible but only over a finite interval Ah of the 
amplitude of the exchange field of the F domains. This width of the temperature induced 
transition increases with decreasing L/C- For cp = 0, the coupling between the exchange 




Figure 5. (Color online) Plot of 1^ versus position for the system with i/^ = 1.0, 
Lf/L = 1/3, h/Tc = 5.0, ip = 0.0, t = 0.1, and QLf^tt; 2tt; Stt and 4tt. Solid, dashed, 
dotted, and dashed-dotted lines, respectively, represent the results for QLf^tt; 27r; Stt 
and 47r. (b) The same as (a) but for ip = n. 



fields of tlie two domains is found to be always antiferromagnetic in an asymmetric 
structure, which is very similar to the symmetric case. 

3.2. S IF 1DWF2S2 junction with Neel domain wall 

The superconducting spin Josephson effect described above may take place in F contacts 
with a more complex profile of the exchange field vector. An interesting case is a finite 
width F domain wall between two domains, where the exchange field has a continuous 
spatial rotation between two homogeneous F domains. Here, we present the results of 
our calculation of the spin supercurrent and the spin transfer torque for a Neel domain 
wall junction, which is shown schematically in Fig.[T]D. Note that the spin transfer torque 
acting on the local magnetization on node i is obtained from r^i = Izi,i+i + Izi,i-i, as we 
have shown earlier. We model the local angle of the exchange field vector with respect 
to the |/-axis to vary as a{x) = {QLp){x / Lp). 

We have obtained the position-dependent spin supercurrent, which is perpendicular 
to the plane of the magnetization and flows in the homogeneously and inhomogeneously 
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magnetized parts of the system. To the best of our knowledge, there have been no 
previous studies investigating the position dependence of the spin current and the spin 
transfer torque. Figures [5](a,b) show the behavior of the spin transfer torque versus 
position for a system with Lp/L = 1/3. We can see that the system has an interesting 
behavior depending strongly on the values of QLp. As we see in Fig. [5K, the spin 
transfer torque penetrating into the homogeneous ferromagnets becomes of negligible 
constant value. This shows that the spin supercurrent in Fi and F2 has a linear 
position dependence when = 0. In particular, this is true while the value of the 
spin supercurrent is comparable to the one in the nonhomogeneous parts. 

Also Fig. [5b shows that the penetrating spin transfer torque in the homogeneous 
parts is nearly zero for QLp = tt, Stt and is much smaller than the one in the domain wall 
region for QLp = 27r,47r when = vr. In addition, we note that the the spin transfer 
torque has always a symmetric position dependence around x = L/2, whereas the 
spin supercurrent always shows an asymmetric behavior. Our calculations demonstrate 
further that the behavior of the spin current and spin transfer torque versus position for 
a system with QLpTc and = vr is similar to that of the system with QLp = (n + l)7r 
and = 0, see Fig. [5l This observation of a symmetry between the magnetic winding 
number and the superconducting phase is presently not fully understood, but will be 
the subject of future research. 

Finally, we have also studied Josephson systems without the homogeneous 
ferromagnetic parts Fi and F2, which corresponds to Lp = L. In this way we would 
like to check how the spin currents in the homogenous part, which are appreciable in 
size, but give rise to a negligible spin-transfer torque, influence the spin-torque on the 
magnetization texture. We show the corresponding dependence of the spin supercurrent 
on the phase difference and the position in Fig. [6] for a wall with rotation angle tt. The 
first thing to note, is that we observe a qualitatively similar phase dependence as in 
the case with two homogeneous ferromagnets. In particular, the Iz{4') relation is always 
a symmetric. Also, for a specific position, the sign of the spin transfer torque can be 
changed by varying the phase difference between the superconductors. In addition, 
for a simple Neel domain wall, when QLp = tt, we see that the behavior of the spin 
transfer torque can be approximately described by the relation cos[{x/L){QLp) + 7i/2] 
for = 0, which turns into cos[3{x/L){QLp) + 7r/2] for = vr. These findings have 
the same symmetry, as we mentioned in the previous discussion for a system with 
homogenous ferromagnet attached. Also, this result shows the strong dependence of 
the spin transfer torque on the phase difference between the superconductors, which 
can be interpreted as a direct interplay between the induced spin-singlet correlations of 
the superconductors and the inhomogeneity of the magnetization in the domain wall. 
In fact this interplay leads to the generation of the triplet correlations in the contact 
region, whose inhomogeneity drives the spin supercurrent and therefore the spin-transfer 
torque. 

Figures [71^a,b) show the spin transfer torque as function of the phase difference, 0, 
and the wave vector, QLp (which is more or less given by the number of 180 degree 
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Figure 6. (Color online) Spin transfer torque Iz versus position and phase difference 
between the superconductors for the Josephson system without any homogeneous 
ferromagnet, Lp/L = 1, while L/^ = 1.0, h/T^ = 5.0, t = 0.5, and QLp=tt. 

domain walls). In Figs. [71^a,b) the behavior for different positions x/L = 1/2 and 
X = L/6 are shown. For the position closer to the superconductor, x/L = 1/6, the spin 
transfer torque goes to zero much faster for larger wave vectors than in the middle of 
the ferromagnet at x = L/2. Hence, we see that the suppression of the spin transfer 
torque strongly depends on the position in the domain wall. Also, Fig. [71^a,b) shows that 
the spin transfer torque oscillations versus wave vector strongly depend on the phase 
difference and also on position. While they start for small QLf in exactly opposite 
fashion for both positions, the oscillations in the middle of the domain wall are well 
behaved also for large QLp, whereas the behavior is more complex in the ferromagnetic 
part close to the superconductor. In that limit no well defined oscillation period can 
be identified. Furthermore, while the direction of the spin supercurrent is fixed, we 
find that the sign and magnitude of the spin supercurrent and spin transfer torque can 
be modulated by changing the wave vector or the phase difference. We have obtained 
that superconductivity in coexistence with non-coUinear magnetism can be used in a 
Josephson nanodevice to create the tunable spin supercurrent which acts as a spin 
transfer torque on the junctions magnetization. 

4. Conclusion 

We have studied the Josephson effect in a diffusive contact consisting of 
two ferromagnetic domains with noncollinear magnetizations which connects two 
conventional superconductors. Using quantum circuit theory, we have shown that 
the spin supercurrent will flow through the contact, due to the generation of an 
inhomogeneous spin-triplet superconducting correlations. The polarization of the spin 
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Figure 7. (Color online) (a) Spin transfer torque of Iz versus QLp and cj) for the 
system with x/L = 1/6,L/^ = 1.0, Lp/L = 1.0, h/Tc = 5.0 and t = 0.1. (b) Same as 
the previous one but for x/L = 1/2. 

supercurrent is directed normal to the plane of two magnetization vectors in a way that 
the resulting spin transfer torque intends to align the magnetization of the two domains. 
The spin and charge current-phase-angle relations obey the specific symmetry relations 
versus the phase difference ip and the misorientation angle 6. Whereas, the charge 
supercurrent satisfies the odd-even relationship on the cp and 6, the spin supercurrent 
is an even-odd function of and 6. From these relations we have predicted a transition 
between and tt Josephson coupling by varying the misorientation angle 6. We found a 
transition of the favorable configuration of the domain, from antiparallel to the parallel 
as the exchange field of the domains increases. This transition occurs for asymmetric 
systems with = tt. Also, the domains in the symmetric systems settle in a parallel 
configuration when (p = tt. 

We have further discussed the generation of the spin supercurrent in magnetic 
contacts with more complex configuration of the direction of the magnetization vector. 
For a domain wall between two domains with antiparallel magnetizations, we have shown 
that the profile and penetration depth of the spin supercurrent are highly dependent 
on the number of the rotations that the magnetization vector has undergone across the 
domain wall. In particular, we show that, while the direction of the spin supercurrent 
is always perpendicular to the plane of the magnetization vectors, the sign and the 
magnitude of the spin supercurrent strongly depends on the phase difference between the 
superconductors and the value of the wave vector. We present a Josephson nanodevice 
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which can be used to create a controllable spin supercurrent and spin transfer torque. 
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